Quantum electrodynamics near anisotropic polarizable materials: 
Casimir-Polder shifts near multi-layers of graphene 
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In a recent paper we have formulated a theory of non-relativistic quantum electrodynamics in 
the presence of an inhomogeneous Huttner-Barnett dielectric. Here we generalize the formalism to 
anisotropic materials and show how it may be modified to include conducting surfaces. We start 
with the derivation of the photon propagator for a slab of material and use it to work out the energy- 
level shift near a medium whose conductivity in the direction parallel to the surface far exceeds that 
in the direction perpendicular to the surface. We investigate the influence of the anisotropy of the 
material's electromagnetic response on the Casimir-Polder shifts, both analytically and numerically, 
and show that it may have a significant impact on the atom-surface interaction, especially in the 
non-retarded regime, i.e. for small atom-surface separations. Our results for the energy shift may 
be used to estimate the Casimir-Polder force acting on quantum objects close to multilayers of 
graphene or graphite. They are particularly important for the case of trapped cold molecules whose 
dispersive interactions with surfaces often fall within the non-retarded regime where the anisotropy 
of the material strongly influences the Casimir-Polder force. We also give a formula for the change 
in the spontaneous decay rate of an excited atom or molecule near an anisotropically conducting 
surface. 

PACS numbers: 31.70.-f, 41.20.Cv, 42.50.Pq 



I. INTRODUCTION 

In a recent paper [l[ we have described a microscopic 
model for the interaction of the quantized electromag- 
netic field with an absorptive dielectric body. The basic 
idea behind the model, first introduced in [2(, is that a 
linearly responding dielectric material can be perceived 
as a continuum of damped quantum harmonic oscillators, 
making up a Huttner-Barnett dielectric. Then the clas- 
sical response of the material to an externally applied 
electric field is described by the Drude-Lorentz permit- 
tivity, 



the reservoir responsible for the absorption, "Hp_R cou- 
ples the polarization field and the reservoir, Hp_em cou- 
ples the polarization and the electromagnetic fields, and 
ris arises from the frequency-shifts due to these two cou- 
plings. In order for Eqs. ©-([7]) to be meaningful the 
following equal-time commutation relations are required: 
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In Ref. [l[ we have demonstrated that the Hamiltonian 
density of such a system can be written as 
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The field D(r) = e E(r) + ,g(r)X(r) is the divergence- 
free displacement field. Its appearance is a consequence 
of the multipolar coupling between the electromagnetic 
field and the polarization field X(r) @. The polarization 
field is in turn coupled to a continuum of bath oscilla- 
tors Y„(r) with variable inertia p v leading to the absorp- 
tion of radiation in the model. The position-dependent 
coupling function g(r) characterizes the volume of space 
where the matter-radiation interaction is switched on. It 
is equal to unity in the volume of space occupied by the 
dielectric and zero otherwise. Therefore it determines the 
shape of the polarizable body in question, cf. Eq. ([lj. 
For a more detailed description of the model we refer the 
reader to Ref. [l|. 
The Hamiltonian 



H = J d 3 r {Hem + n P +n R 



H 



P-R 



H 



P-EM 



\ r rf ^X 2 (r) + i^X 2 (r), (7) 
^ Jo 1 e o 



where Hem describes the dynamics of the electromag- 
netic field, Hp that of the polarization field, Hp that of 



can be in principle diagonalized, see e.g. j4[, and the ex- 
act electromagnetic field operators can be written down 
for some simple geometries of the dielectric body. How- 
ever, as we have shown in [l[, it is also possible to sidestep 
the somewhat complicated diagonalization procedure and 
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instead use the well-known methods of quantum field the- 
ory to directly obtain the exact Feynman propagator for 
the displacement held. The set of Eqs. (O-® define an 
interacting quantum field theory and a rather simple one 
where all the couplings are bilineax. The chief quantity 
of interest is the propagator of the displacement field, 
which couples to any additional electric dipole moments 
in the system and as such can be understood as medi- 
ating interactions between atoms and molecules through 
their (fluctuating or permanent) electric dipole moments. 
With knowledge of this propagator one can calculate 
various experimentally measurable quantities, e.g. the 
Casimir-Polder shift in a neutral atom or molecule, in 
much the same way as this is done in free-space quantum 
electrodynamics. The Feynman propagator is an auxil- 
iary mathematical construct which happens to appear in 
the diagrammatic formulation of pcrturbativc quantum 
field theory and in general it does not necessarily have a 
clear physical interpretation However, in our model 
we may still associate the propagation of photons with 
the displacement field propagator defined as 

D y (r,r;t,t / ) = -|(0|T[A(r,t)D J -(r / ,O]|0) , (11) 

and therefore we shall call (r, r; t, tf) the photon prop- 
agator. The state |0) in Eq. (fTTj) denotes the exact 
ground state of the interacting system, and Di(r,t) is 
the displacement field operator in the Heiscnberg picture 
whose dynamics is governed by the full complement of 
Hamiltonian densities ©-([7]), that is to say, including 
all the couplings. Using basic diagrammatic techniques 
we have shown in Ref. [l[ that the Dyson equation for 
the photon propagator can be written as 

D«(r,r'; W ) = 4 0) (r-r'; W ) 
+ fd^gir^Dlfir-^D^r'^) (12) 

where D^\r— t';uj) is the free-space photon propagator 
and K(tu) is a frequency-dependent function related to 
the permittivity of the medium (see Sec. IV of Ref. [1| 
for details) . The benefit of working with propagators is 
that with Feynman diagrams one has at one's disposal 
a very efficient bookkeeping device for keeping track of 
higher-order contributions to perturbation series, which 
otherwise is a non-trivial task. As an example of what 
quantum field theory can do for atomic physics consider 
the van der Waals interaction of two identical atoms with 
one of them excited. If the interaction energy is calcu- 
lated by using fourth-order time-independent perturba- 
tion theory one needs to consider twelve graphs and the 
same number of contributing terms , but in quantum 
field theory the calculation of the same quantity requires 
only two Feynman diagrams Q. 

It is well known that there are materials whose re- 
sponse to electromagnetic fields is anisotropic i.e. which 
depends on the direction of the applied field. Imagine 
e.g. a slab of the material, cf. Fig. [TJ which in the direc- 
tion parallel to its surfaces behaves as a good conductor 
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FIG. 1. (Color online) The Feynman propagator for the dis- 
placement field, given in Eq. (|48[) . describes the propagation 
of photons near a slab of absorptive material with anisotropic 
and frequency-dependent permittivity e a (to). It is calculated 
exactly starting from the quantum model for the light-matter 
interactions described in Sees. HlfcHTl 

but in the direction perpendicular to its surfaces is a di- 
electric (or weak conductor) Q. An example of such a 
material would be a single layer or multi-layer films of 
graphene, which have recently started to play a signif- 
icant role in nanotechnology, nanoelectronics and even 
experimental cold-atom physics. Thus, it seems to be of 
great importance to have appropriate tools at hand for 
the calculation of radiative corrections to properties of 
quantum systems, like atoms and molecules, near such 
anisotropic interfaces. The main goal of this paper is to 
provide those tools. In what follows we are going to show 
how the Huttner-Barnett model described above can be 
generalized to anisotropic dielectric media. As an exam- 
ple we will calculate the Casimir-Polder force acting on 
an atom close to an anisotropic dielectric slab. This pa- 
per will heavily rely on previous work reported in Ref. [l| , 
from where we will be quote results instead of re-deriving 
them. So, this paper should be read in conjunction with 
Ref. [l| , in particular with Sections III and IV thereof. 

II. PHOTON PROPAGATOR NEAR 
ANISOTROPIC MEDIA 

For an anisotropic dielectric slab, as depicted in Fig. 
[TJ the permittivity needs to be generalized from a scalar 
to a tensor 

/ e\\(z,u) \ 

e(z,u) = e\\(z,(j) , (13) 
\ e±(z,u) J 

where ej_ and ey are still of the Drude-Lorentz form of 
Eq. ([JJ and, for a slab, depend only on z coordinate. In 
order to achieve a dielectric response of the form (|13[) in 
our model we introduce an anisotropic polarization field 
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and replace Eq. ([3]) with 

Pff(r) 1 



Hp = 



2M\ 



X||w4||X|(r 



+ ^ + \ M ^ Xl{T) - (14) 

To account for the slab geometry we fix the coupling 
function in Eqs. ([6]) and (|Tj) to 



9(r) = 1 



(15) 



as appropriate for a dielectric slab of width L centred on 
the z = plane, as shown in Fig. Q] Here 9(z) is the 
Heaviside step function. 

In order to derive a generalization of the integral equa- 
tion (|T2"j) that reflects the anisotropy of the dielectric, one 
needs to successively integrate out the matter degrees of 
freedom from the Hamiltonian density. This can be done 
along exactly the same lines as described in Sees. Ill & 
IV of Ref. [l[ . The only difference is a change in the po- 
larization field propagator on account of the Hamiltonian 
density having changed from Eq. ([3]) to Eq. f| 14[) . We find 
that the dressed polarization field propagator is still di- 
agonal, but its diagonal elements arc no longer equal as 
they were for an isotropic medium. Now we have 



with 



eo 



K {1 (uj) 
A",, (ui) 
K±(u) 



& (w) = 1 + 



$< 3 >(i 



(16) 



2 

PfT 



(17) 

where a = {||, _L}. Following the same steps as described 
in Sees. Ill & IV of Ref. [lj, we derive the generalized 
integral equation for the photon propagator, 

AKr,r» = Df (r-r';w) 

+ \ /dVsCri^Cr- n;o;)X 3 - J fc(w)D w (ri > r';w) (18) 
e o J 

with the free-space photon propagator given by 



r :u 



(2 



27r) 3 J lo 2 - q 2 + it] 



(19) 

The solution of Eq. (|18p . with g(r) as in Eq. (fTS"]) , is going 
to yield the photon propagator for the displacement field 
near an anisotropic dielectric slab, and with some skill 
it can be obtained by iteration. We are going to rely 
on a similar trick to the one used in Ref. [l| for solving 
the equivalent integral equation for the geometry of a 
dielectric half-space. Equation (fT5| . as it stands, cannot 
be iterated directly but can be reformulated in such a 



way that iteration is possible. For the slab geometry the 
matter-radiation interaction Hamiltonian reads 



^P-EM — 



1-0 -z-- - 



L 

z 

2 



X(r)-D(r). 

(20) 



There are several ways of splitting this Hamiltonian into 
an unperturbed part, for which the photon propagator 
is known, and one to treat as a perturbation, and dif- 
ferent choices lead to different integral equations but for 
the same photon propagator. Equation (|18p has followed 
from taking the whole of Eq. (|20p as the perturbation 
and the unperturbed part of the Hamiltonian as that of 
the free electromagnetic field. However, we are free to 
split up the Hamiltonian differently, provided we know 
the photon propagator for the unperturbed part of the 
Hamiltonian. Let us do this in two different ways. First 
we split the interaction Hamiltonian so as to include the 
polarization field for a dielectric half-space to the right 
of —L/2 in the unperturbed Hamiltonian, 



EM 



1 



1-0 -z 



H 



P-EM 



1 

— ( 



L 



--]X(r)-D(r). 



X(r)-D(r), 
(21) 



With this choice of splitting, the integral equation for the 
photon propagator reads 

A/(r,r'; W ) = ^ +) (r,r'; W ) 

(d?Y X e \z! - ^ D<+>(r, n;w)lT ifc (w)D H (r 1> r / ;w) 

(22) 

where (r, r'; lj) is the photon propagator for the case 
of an anisotropic dielectric half-space occupying the re- 
gion z > —L/2. Next we reverse this choice and split the 
interaction Hamiltonian so as to include the polarization 
field for a dielectric half-space to the left of L/2 in the 
unperturbed Hamiltonian, 



H 



(o) _ 



1 

eo 



1 - 6 




H)] 







X(r) • D(r), 

(23) 



With this choice of splitting, the integral equation for the 
photon propagator reads 

D a (T,r';oj) = D^\r,v';uj) 



z i - \ ) D l ) (r,r 1 ;o;)X J - J fc(w)D w (r 1> r';w) 



(24) 



where (r, r'; uj) is the photon propagator for the case 
of an anisotropic dielectric half-space occupying the re- 
gion z < L/2. 
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To keep the notation concise, we will from here on work 
with quantities that have been Fourier transformed in the 
direction parallel to the surface, e.g. with 

Dyfoa^qil.w) = /"d 2 R||e- iq ll- R llAj(R||,2,^;a;). 

(25) 

where Ry = — r'| . For notational clarity we will also 
suppress the dependence of the propagators on the fre- 
quency u> and on the component qy of the wave vector 
parallel to the interface. 

We proceed by substituting Eq. ([24]) into the right- 
hand side of Eq. (f2"2"j) , which gives 



D ij (z,z')=D% ) (z,z') 



1 



dziD%\z,zi)K 1d {u>)D\J>{z u z!) 

e JL/2 
^ poo p—L/2 

+ — / <fci / tfe 3 -D|fc" ) (^«i)-K'fcl( w ) 
xD { ~J ( Zl , z 2 )K mn {Lu)D nj (z 2 , z'). (26) 



>(-)/ 



Integral equation (f2"6")) is equivalent to Eqs. (|18p. (|22[) . and 
(|24|l in the sense that their solution is the same photon 
propagator for the anisotropic dielectric slab. The crucial 
difference is that the integral equation (|2"6"|) can be solved 
by iteration, whereas the others could not. To show that 
this is so, we need to know the photon propagator near 
an anisotropic dielectric half-space. The latter can be 
calculated by building on the expertise gained from the 
calculation of the photon propagator near an isotropic 
dielectric half-space and reported in Ref. [l|. 



A. Photon propagator near an anisotropic 
dielectric half-space 

Let us consider a half-space that occupies the re- 
gion z < 0. We seek the solution of Eq. (fT8|) with 
g(r) = 6(—z), which can be found by iteration. We only 
outline the calculation in this section, as we are closely 
following the steps of Sec. IV B of Ref. and we shall 
concentrate on the modifications needed to account for 
the anisotropy of the material. The calculations are very 
similar, although a little bit more tedious. For a half- 
space filling the region z < 0, the photon propagator 



satisfies the following set of integral equations 
D u (z,z') = DP(z-z') 



+ 4/ d Zl D^\z- Zl )K lk {^)D kl {z u z'). (27) 



ropa 



D u (z,z')=D%\z-z') 
1 f°° 

—5 / dzxDfhz-z^K^Duiz^z'), (28) 



c o Jo 



were Dw\z — z') and D^\z — z') are the propagators 
in free space and in a bulk (now anisotropic) medium, 



respectively. Equations (|27P and (|28|) combine to yield 



D in (z,z')=D%(z-z') 



- 2 / d Zl Dlf(z - z^K^D^iz, - z') 



L JO 



(0), 



c o Jo 



dz x 



dz 2 D\f \z - Z!)K jk {u) 



xD$\z! - z 2 )K lm {ui)D mn (z 2 ,z') . (29) 

In order to be able to solve Eq. ([2"9"| . we need to de- 
termine the displacement field propagator in the bulk 
anisotropic medium, that is, the quantity D^\z — z'). 
To this end we write down Eq. (fT5|) with g{r{) = 1 and 
Fourier transformed with respect to r — r', 



<%q 2 - Qilj Kjijoj) 
u! 2 — q 2 + ir\ eo 



e — o 



uj* — q^ + vr\ 
(30) 

where we have used the spectral representation of the 
free-space propagator from Eq. (fTUj) . Thus, finding 
is just a question of finding the inverse of the matrix in 
the square brackets on the left-hand side. This matrix 
may be written as 



17] 



e, A (qK(q)- 



(31) 



where we have used the completeness property of the 
transverse polarization vectors defined as 

1 

hi 

-(q x k z ,q y k z ,-q L j). (32) 



e TE (q||) = Tr-rfej-^O), 



TM, 



qi , k z 



|qiil w 



The products of the polarization vectors are projection 
operators onto the two orthogonal polarizations, so that 
their repeated application always gives the same polar- 
ization again. We have in particular 

e A (q) ej A (q)A^^)e7(q) e Z(q) = **7V(q)e£(q) (33) 
with 

/ TE ^iiM, 

1 q 2 
With that we are able to find the inverse of the matrix 
(|3"Tj) by writing it as sum over a geometrical series, ac- 
cording to 



(1-0)- 



l + O + O 2 



Equation (|33|l allows us to carry out the infinite summa- 
tion, and we find 



^•(r-r'jw) 



(2tt) 3 

iq- (r— r') 



xe 



yd 3 qQ A (q, W ) e f(q) ej A (q) 
(34) 
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with 



has 



Q TE (q,-) 
Q™(q,w) 



£y (w)w 2 — q 2 ' 



This result reproduces the bulk-medium propagator in an 
isotropic medium if we take = The trans- 

verse electric contribution to the propagator depends 
only on which is not surprising because for this 

polarization the electric field in the direction perpendic- 
ular to the surface vanishes. Therefore, the TE part of 
the propagator is insensitive to the response of the ma- 
terial in the _L direction. The TM contributions to the 
propagator of course depend on both £||(w) and 

The quantity that enters the Eq. (f29|) is the propaga- 



tor ()34p but Fourier transformed back from q z to z — z' . 
Carrying out the q z integration, one finds away from the 
point z = z' 



A\2 



em d )^ zd )e^-^ z > z> 



(35) 



with q* = (qn , k* d ) and 



h™ - 



where the square roots are taken such that their imag- 
inary part is always positive. The medium affects the 
polarization vector of the TM mode which is given by 



,TM 



fall 



> K zd ) 



|qJ M l 



{lxk zd ,q y k zd , — qn) 



An equivalent representation for free-space propagator 
Dij\z — z') will also be needed in the following; it may 
be obtained simply by taking the limit £o-(w) — > 1 in 
Eq. (|35|) . We emphasize that formula (|35j) is valid only 
away from the point z = z', where additional singular 
terms would contribute. We also note that for the evalu- 
ation of the z integrals in Eqs. (|27|) - ((29|) the polarization 
vectors need to be written in terms of differential opera- 
tors acting on e ^ll-(*ll-*ii)+«£il«-*'l. 

Using the bulk and free-space propagators derived 
above, one can verify that for z < and z' > one 



c "'O 



dzi 



dz 2 Dff (z - Zi)K Jk (u) 

O 

-,(0) 



xD^(z 1 -z 2 )K lm (oj)D^ n (z 2 -z') 
= £A D SJ Z -A ( 36 ) 

A A 

This features the Fresnel reflection coefficients at an 
anisotropic dielectric half-space, 



„TE 



TM 



TM 



J.TE ' 
K zd 



(37) 



with k 7 



'■ — q| + ii], which is the z component 

of the wave vector in vacuum. This has a small posi- 
tive imaginary part which originates from the prescrip- 
tion of handling the poles in the free-space propagator in 
Eq. (|19p and serves to render the Fourier integrals mean- 
ingful. 

Repeated use of Eq. ([55]) facilitates the iteration of the 
integral equation (|29p . yielding the expansion 



D\ :l i(z,z') = 



D%{z-z') 



~ I dz 1 D%{z-z 1 )K jk ^)D^ H {z 1 -z') 
e o Jo 



,(0) 



f - ri 



l-r\ 



(38) 



This is the photon propagator for the case z < and 
z' > 0. Substituting the bulk and free-space propagators 
and summing the geometric series gives 



Dij(z, z') 



xe*(q h -k z y k * z '- lk ^ z 



with 



.TE 



jTM 



2y^J)k z 



(39) 



(40) 



The propagator for the case z, z' > can now be derived 
by applying the integral equation (|27[) . It turns out to 
have the same form as that for the isotropic dielectric, 
derived in Ref. [l| , but with the Fresnel reflection coeffi- 
cients now those of Eq. (|3"T|) for an anisotropic half-space, 

D ij (z,z')=D\f(z-z') 



~2kT 



^r A e^(q||,fc z )e^(q||,-fc z ) 



Xf 



ik z (z-\-z') 



(41) 



The photon propagator in coordinate space is obtained 
by inverse Fourier transform, i.e. by inverting Eq. (|25p . 
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B. Photon propagator near an anisotropic 
dielectric slab 



With the photon propagator for a dielectric half-space 
determined, we are in position to work with integral equa- 
tion (|26[) and obtain the photon propagator for the slab 
geometry. Aiming to solve Eq. ([2"B"]) by iteration, we eval- 
uate 



1 

e JL/2 



-L/2 



dz 2 D^\z,z 1 )K kl {u)D ( J(z 1 ,z 2 ) 



;K mn (u;)D%\z 2 ,z') = ^D^(z,z') ( 



A ifc A ,L 

r e zd 



(42) 



where we have restricted z < —L/2 and z' > L/2. To 
arrive at Eq. (|4"2"]l we have used the result ([55)1 . trans- 
lated by L/2 for D^~\z, z 1 ), and reflected and translated 
by -L/2 for D (+ \z,z'). Equation (02]) shows that the 
action of the integral operator in Eq. (I26| amounts to a 



simple multiplication by a factor of ^r A e lfe =d L j . There- 
fore the iteration of Eq. ([26)) yields a geometric series 
that can be summed up to all orders, 



Dx,ij(z,z') 



D 



(+) 



1 f°° 

(z,z')--2 / dziDgifozi) 

e JL/2 



x /u ; ( W )^i i ;.(z 1 ,z') 



i 



(43) 



This is the photon propagator in the slab geometry for 
the case z < —L/2 and z' > L/2, i.e. it describes the 
transmission of photons across the slab. The Z\ integral 
can be evaluated, which allows us to write the final result 
in a more familiar form, 



Dij(z,z') 



1€qUJ 

~^2k7 



Y,T x e x M\\-kz)e]{^,-k z ) 



xe 



ik z (z — z r ) 



(44) 



where we have introduced the Fresnel transmission coef- 
ficients at an anisotropic slab, 



1 r A c i(kj,-k z )L 



1 - r\e 2ik 



(45) 



and r\ are the same single-interface reflection coefficients 
as in Eq. (|37| but we have written the polarization as a 
lower index for notational convenience. In order to obtain 
the propagator describing the reflection of photons from 
the slab, i.e. for z, z' > L/2, we use the integral equation 
([M]) . Substituting Eq. (jH| into the right-hand side of 
Eq. JM]) and using the results ([51]) and JIT]), we nn d f° r 
the propagator in the region z,z' > L/2 

D ij (z,z') = Dl"\z,z') 



1€qCJ 



(46) 

where Df^\z,z') is the (partially Fourier transformed) 



free-space photon propagator of Eq. (|19|) . and the Fresnel 
reflection coefficients i? A at an anisotropic slab are given 
by 



-ik z L 



i 



2„2ife A 



rte 



(47) 



As one might have expected, the photon propagator near 
an anisotropic dielectric slab has exactly the same form as 
the one near an anisotropic dielectric half-space, Eq. (|41|) . 
except for the different Fresnel coefficients which account 
for each particular geometry of the dielectric. 

Eqs. ([44]) and ([46]) form a major result of this pa- 
per; they give the photon propagator for the quantized 
electromagnetic field near a slab of anisotropic dielectric 
material that is capable of absorbing radiation. Note 
that we have derived the propagator starting from a well 
defined microscopic model that is explicitly described by 
a Hamiltonian. We emphasize in particular that we have 
not made any reference to statistical physics for relating 
quantum propagators to classical Green's functions of the 
wave equation as is often done in the literature. Although 
the propagators we have derived describe the propagation 
of photons through inhomogeneous dielectric media, we 
did not have to explicitly invoke Maxwell's equations and 
continuity conditions following from them in the course of 
the derivation because this information is automatically 
included in the dynamics prescribed by the Halmiltoni- 
ans. Also, note that throughout the derivation we have 
dealt only with physical electromagnetic fields and not 
with potentials that would have rendered our derivation 
gauge-dependent. For future reference we provide the 
complete result in coordinate space. For the source point 
to the right of the slab, z' > L/2, the photon propagator 
outside the slab reads 



Oij(r,r';w) = 



rf 2 qiie 



(2ny 



lc iir( r ii- r !i). 



2A- 



^TV(q|h-^)e, A (q||,-fc z )e* 



for z < — - 



(2ny 



iqn-( r ii- r ii) 



u> 

2k; 



Y / R X e X (q h k z )e x (q ll ,-k z )e lk ^ +z ' , > for z > 



(48) 
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III. ATOM-DIELECTRIC INTERACTION AND 
ENERGY-LEVEL SHIFTS 

To study the interaction between an atom or molecule 
with an anisotropic dielectric slab, we choose a one- 
electron atom, either in its ground or in an excited state, 
interacting with a quantized electromagnetic field whose 
dynamics is governed by Eqs. ([2|h([Z|)- We treat the 
atom using Schrodingcr quantum mechanics but employ 
methods of second quantization. With the atomic en- 
ergy levels known, the second-quantized Hamiltonian of 
the atom Ha can be written as 

H A = Y J E n c\ l c n , (49) 

n 

where the creation and annihilation operators c\ and Cj 
create an excitation in the atomic energy level i and de- 
stroy one in j, respectively, and satisfy standard anti- 
commutation relations {ci,cj} = 6{j. The atom is cou- 
pled to the field by a multipolar Hamiltonian, and the 
leading contribution to the interaction energy between 
atom and dielectric comes from the interaction between 
atom's electric dipole moment and the electromagnetic 
displacement field, which mediates the atom-surface in- 
teraction. In its second-quantized form the interaction 
Hamiltonian reads 

Ha-em = -- £ 4 Cj (i\p\j) • D(R), (50) 
y 

where —e(i\p\j) are the matrix elements of atom's elec- 
tric dipole moment and R is the position of the atom, 
which we always assume to be at least few Bohr radii 
away from the surface so as to be allowed to neglect di- 
rect wave function overlap between the atom and the di- 
electric medium and surface chemistry. The energy-level 
shift in the atom arises as radiative correction to the 
self-energy of the atomic electron and may be calculated 
using the standard tools of quantum field theory, if one 
exploits the fact that the poles of a Feynman propagator 
correspond to the energy spectrum of the excitations [t| . 
Here we are after the atomic energy-levels which give rise 
to poles in the Fourier transform of the time-dependent 
propagator that describes the state of the atomic electron 
and in our formalism is given by 

Glrn{t,t') = - l -(n\T[ci{t)cl{t')} |fi). (51) 

Here Ci(t) is a time-dependent operator in the Heisen- 
berg picture, |f2) is the exact ground state of the system, 
and T is the time-ordering operator. The coupling (f5"0"l) 
between the atom and the dressed electromagnetic field 
changes the analytical structure of the atomic propagator 
(|5ip in the interacting as compared to the non-interacting 
case and shifts its poles and thereby the energy-levels 
of the atom. One part of this shift is the well-known 
Lamb shift, which is the same everywhere in space, and 
the other part depends on the distance of the atom to 



the surface because the interaction Hamiltonian Ha-em, 
Eq. (|50p . contains the displacement field operator which 
takes into account the presence of the dielectric via its 
dynamics as prescribed by the Hamiltonians in Eqs. @- 
([7]) and is therefore governed by the photon propagators 
in Eqs. (|4Tj) or (|48|) for the half-space or the slab, respec- 
tively. 

The formulae for the energy-level shift of an atom in- 
teracting with the electromagnetic field as a functional of 
the photon propagator has been derived step by step in 
Ref. [l| . The energy- level shift splits into two distinctive 
parts, 

A£[ on = AS, + AS* , (52) 

where the superscript 'ren' indicates that we are inter- 
ested only in the rcnormalizcd energy-level shifts and 
have subtracted the free-space Lamb shift. In practice 
this renormalization procedure is very easily achieved 
thanks to the fact that the photon propagator out- 
side the dielectric always splits into the free-space part 
D*j? '(r — r';uj) and the part that describes the reflection 

from the surface D$(r, r'; u>). Therefore, we just need 
to disregard the free-space part of the photon propaga- 
tor and work only with the reflected part. With this in 
mind, the energy shift ([52)1 can be written as [l[ 

A£=-^5>£j a / ^-^VZ^(R,R;zO,(53) 




where l/ij^il = K m lM fe |*)l are the matrix elements of the 
k-th component of the electric dipole moment operator 
and oj m i = ui m —uii is the frequency difference between the 
unperturbed atomic energy levels \m) and \i). The sums 
in Eqs. (j53|) and (|54|) run over the Cartesian components 
of the electric dipole moment operator k = {x, y, z} and 
over all atomic states |m) excluding the state \i) that we 
are calculating the energy shift of. The sum over m is in 
practice limited to one or a few states to which there are 
strong dipole transitions from the initial state The 
shift ASi affects all atomic states but the contribution 
AS* arises only if \i) is an excited state. Equations analo- 
gous to Eqs. (|53| and (f54|) can also be derived by different 
methods e.g. by linear response theory (ToL or by a 
phenomenological noise-current approach to macroscopic 
quantum electrodynamics [T2|. The shift ASi is always 
real because the photon propagator is real at complex 
frequencies. However, AS* is complex and contains cor- 
rections to the spontaneous decay rates of excited states. 
One has 

AEi = Re (A£,[ on ) 

ATi = ~Im(A£*) (55) 

where AEi are the energy-level shifts and Ar^ are the 
changes in decay rates. 
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The photon propagators that we have derived in Scc.HIl 
are formally of the same form as the photon propagator 
that was derived in Ref. [l| for the case of an isotropic 
dielectric half-space. The difference lies only in reflection 
coefficients. For this reason most of the formulae for the 
atomic energy-level shift derived in Ref. [l| can be in- 
stantly generalized to the case of the anisotropic dielec- 
tric media studied here. First we are going to study the 
interaction between a neutral atom and an anisotropic 
dielectric half-space. Because the half-space geometry 
is relatively simple, we are able to study the effect of 
the anisotropy on the Casimir-Polder interaction in some 
detail. Then we go on to explore the physically more 
relevant case of a medium of finite thickness. But as the 
case of an anisotropic dielectric slab results in more com- 
plicated formulae for the energy-level shift, we will give 
fewer analytical and more numerical results than for the 
half-space. 



IV. CASIMIR-POLDER SHIFTS NEAR AN 
ANISOTROPIC HALF-SPACE 

In this section we are going to derive a formula for the 
energy level shift of an atomic electron due to the pres- 
ence of an anisotropic polarizable half-space. For this we 
focus on ground state shifts, but in Appendices lAl and IBl 
we also provide expressions for the shifts of excited en- 
ergy levels and for the changes in the spontaneous decay 
rates due to the interaction with the half-space. 



A. Ground state shifts 

In order to obtain the general expression for the en- 
ergy shift, we substitute the reflected part of the photon 
propagator (|4Tj) into Eq. (|53p and obtain 



AE„ 



8tt 2 



-y 



dkk 



-2%/W+wZZ 



+ UJ^g 



y/k 2 + w 2 



{ [(k 2 + w 2 )f™ - w 2 -TEj |/4 y 2 + 2fc¥ M |M i j2 j ; 

(56) 



\t4ni\ 2 



| M LI 2 - There- 



with the abbreviation l/i^J 2 

flection coefficients are defined in Eq. (|37|) : in terms of 
the new variables used in Eq. (|56|) they read 



-TE 



-TM 



Vu 2 + k 2 - yJt\\L0 2 + k 2 

\/ui 2 + k 2 + ^e^ui 2 + k 2 ' 



v /ejjeXVw 2 + k 2 - \/e ± uj 2 + k 2 
^elv'uj 2 + k 2 + y/e±uj' 2 + k' 2 ' 



(57) 



(58) 



where all dielectric functions are evaluated at imaginary 
frequencies, that is, e a = e a (iuj). Note that compared 
to Eq. (|37[) we have replaced the previously unnamed 
function ^(w) with the dielectric constant e a (u)) because 



both functions coincide when evaluated at imaginary fre- 
quencies, cf. Eq. p7|) . 

The expression for the energy shift in Eq. ([55)1 is not 
necessarily the most convenient for asymptotic or numer- 
ical analysis. To find an alternative expression, we go to 
polar coordinates, uj = xoj mg cos <p, k = xui mg sin0, and 
then set cos0 = y. After short calculation wc obtain 



1 



87r 2 e Z 4 



E: 



J mg 



J-I..-L |2 
mgl 



with \/j)L g \ 2 



/x^,J 2 and 



F" 



\^mg\ 

1 (■?. . \4 



(Zoj mg ) 



dxx 3 I dy 
o Jo 1 + £ V 



TM -2a 



(59) 



-,Zx 



(60) 



,Zx 



(61) 



In these variables the reflection coefficients read 



1 + 1 



1+ sjy 2 {e\\ -1J + T 



jTM = 



V^T- y/y 2 (e± - 1) + 1 
V^T+v / y 2 ( e i- 1 ) + 1, 
where e a = e a {ixyu mg ) 1 that is, for example, 



(62) 



e\\(ixyuj mg ) =l + - 



(ujpn/uj mg )' 2 



x 2 y 2 + (cj T ||/u; mg ) 2 + 2xy(^\\/uj mg ) ' 



Thus all frequencies entering Eq. (|59p arc measured in 
units of the atomic transition frequency uj mg and the 
argument of the exponential Zuj mg — 2TrZ/X mg is the 
atom-surface distance measured in units of atomic tran- 
sition's wavelength A mg . Note that the sum over m in 
Eq. (|59p also contains the dipole matrix elements, whence 
in practice it is dominated by just one or a few states with 
strong dipole transitions from the initial (ground) state 
\g). We shall refer to u> mg as the frequency of a typical 
atomic transition. 

The expression for the energy shift given in Eq. (|59p (or 
equivalently in Eq. (|56p ) is general in the sense that en (cj) 
and ej_(w) are arbitrary as long as they can be repro- 
duced by the oscillator model that we have started with. 
Note in particular that this includes dielectrics with more 
than one absorption line, as discussed in Refs. [H, 0- In 
the following we would like to consider the asymptotic 
behaviour of Eqs. (|56p and (|59p in a few physically im- 
portant limits. 



1. Nonretarded limit 

First we consider the non-retarded (or electrostatic) 
limit where electromagnetic interactions are instanta- 
neous and the speed of light is infinite, c — > oo. Having 
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worked hitherto in natural units where c = 1, we recall 
that for the frequency ui to have the same dimensions 
as the wave vector k it needs to be multiplied by 1/c. 
Restoring the missing factors of 1/c in the appropriate 
places and taking the limit c — > oo, we can then perform 
the k integration in Eq. ([SS]) . and the final result can be 
cast in the form 



a pnonrct ^ \ 1 

9 ~ 327r 2 e 2 3 ^ 



rng I 



2|m^I 2 



duo 



'rng (^)e_L (^) ~ 1 

w2 + w m S v /e||(iw)e_ L (iij) + l' 



(63) 



Thus in the nonretarded regime, the ground-state shift 
caused by the anisotropic medium depends on the dis- 
tance Z of the atom from the surface as Z~ 3 (as it also 
would for an isotropic dielectric half-space or a perfectly 
reflecting surface because in all these cases the underlying 
interaction is an electrostatic dipole-dipolc interaction in 
the same geometry). The coefficient of Z~ z arises as an 
integral over the atomic polarizability along the imagi- 
nary frequency axis with a factor 



,J e\\(iui)e±(iuj) — 1 
y/e\\ (iu)e±(iu) + 1 



(64) 



which is a generalization of the electrostatic image fac- 
tor that would have arisen if one were to determine 
the Green's function of the Poisson equation for an 
anisotropic dielectric half-space. We note that setting 
en = ej_ and neglecting the damping leads to the result 
reported in Ref. [lH where the shifts were calculated by 
explicitly quantizing the Maxwell field coupled to a half- 
space filled with a plasma. Formula (f6"5")) is also valid for 
dielectrics or, in fact, for any material whose dielectric 
function can be described by the Lorentz-Drude oscilla- 
tor model. In other words, the electrostatic limit c — > oo 
does not interfere with the limit cjt — > which corre- 
sponds to the dielectric —> conductor limit. This is not 
always the case and in general the order in which various 
limits are taken matters and great care must be taken in 
asymptotic expansions of energy-level shifts [ijj]. 

The limit c — > oo corresponds to the atom being very 
close to the surface, by which we mean that Zuj mg <C 1 
where u> mg is the frequency of the dominant dipole tran- 
sition in the atom. The parameter 2Zui mg /c (if we re- 
store the factor of 1/c) serves as a measure of how fast 
the atom evolves as compared to the time taken by a 
photon for one round trip between the atom and the sur- 
face. If Zu> mg <C 1, the atom does not evolve appreciably 
while the photon travels to the surface and back, which 
is equivalent to taking the speed of light as infinite. 



2. Retarded limit 

In the opposite regime, if Zuj mg 3> 1, retardation is 
crucially important and well-known to have a significant 



impact on atom-surface interactions [15(. The retarded 
limit Zuj m g — > oo is not interchangeable with either the 
dielectric — > conductor limit, wt —> 0, nor with the limit 
of no damping, 7 — > 0. Therefore, when working out the 
asymptotic expansion of the energy shift in the retarded 
limit, we need to specify a definite form of the electro- 
magnetic response of the material from the outset. Let 
us first examine the case of an anisotropic conductor with 
losses, where wc take 



e ± (u) = 1 + - 



e,,H = 1 + 



ui 27j_ — iui 
to 27|| — iuj 



^i + z^M 



= l + i 



(65) 



where a is the usual Drude conductivity. For Zuj mg 3> 
1 the x integral in Eq. ([59")) is strongly damped by the 
exponential. Then, according to Watson's lemma fl6| . 
we obtain an asymptotic expansion of the integral by 
keeping the exponential, expanding the remaining part 
of the integrand in a power series around x = + , and 
carrying out the x and y integrations analytically. We 
find that the asymptotic expansion of the F p coefficients 
entering Eq. (|59[) can be written as 



F" 



9/2 



y/(T\\(0)Z <T\\(0)Z 



(66) 



with <Tii(0) = wpii/27 and 



C 9/ 2 ~16V2' G9 / 2 "T2V2' 65 ~4' ° 5 ~2- 



The first term in Eq. (|66|) reproduces the result of Casimir 
and Polder for a perfectly reflecting mirror [l5l ]. This is 
what one would expect because in the retarded regime it 
is the static response of both the atom and the material 
that matters the most fiol ] . In our model of an anisotropic 
conductor we have perfect screening in the static limit, 
i.e. e(uj) — > 00 for lo — > 0, just like for an isotropic ma- 
terial. That is why the leading terms of the asymptotic 
expansion (|66|) do not depend on a± and are exactly the 
same as for an isotropic medium with mi = <r±. In other 
words, to leading-order, the anisotropy of the conductor 
does not matter in the retarded regime. 

It is worth noting that the expansion (|66[) is not com- 
patible with the vacuum limit wpii — > 0. In this limit 
we would expect the shifts to vanish, but this is clearly 
not the case for formula (|66|) as it contains the perfect- 
reflector terms which are independent of mi (0). Not only 
that; the subsequent terms in the expansion even diverge 
in the vacuum limit. This is because we have no right 
to expect the retarded limit Zuj mg 3> 1 to be compatible 
with the low-conductivity limit Zu)p\\ -C 1, as has been 
noted before [l3[. 

Note also that in the limit 711 — > of no damping only 
the leading-order term survives. One might be tempted 
to think that in the absence of damping any corrections 
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to the perfect-reflector behaviour of the energy shift are 
of higher order than (Zuj mg )~ 5 , but this is not correct. 
The no-damping limit 711 — s- is not interchangeable with 
the retarded limit Zio mg — > 00, and this becomes more 



with 



apparent if we derive the next term of the expansion (|66p , 
which turns out to diverge in the limit of no damping. 
Therefore, the undamped case 7|| — > must be treated 
separately by first setting 7|| — in Eq. (|65p and then 
repeating the steps leading to Eq. (|66|) . Then we get for 
an undamped conductor 



(67) 



with Cir = 4/5 and C\ = 2. 

Since the asymptotic expansion (|66|) applies only to 
anisotropic conductors, i.e. only to materials whose elec- 
tromagnetic response is described by Eq. (|6"5"|) . let us con- 
sider for comparison a case where the material behaves as 
a non-dispersive dielectric in the direction perpendicular 
to the interface and as a lossy conductor in the direction 
parallel to the surface. In other words, we take 




eii(w) = 1 + i 



(68) 



Then, repeating the steps described in the paragraph be- 
low Eq. (|^5|) . we find that in the retarded regime the 
coefficients F p entering Eq. (|59p are given by 



F" 



a 



9/2 



c 



y^To)z <7||(o)z 



, (69) 



with 



a 



9/2 



°9/2 



21 35 1 

64 + 64^7 2 1 



7T 

7T7 



2 8n 2 ± 



1 3 7 

2'4'4' 



1 



3n± 



2F1 I 



1 

1 3 

2'4 



— 9F-1 

2 



C = - 



3 3n^ 



2 



4n 2 ± 



As one can see, the coefficients Cg',^ are rather compli- 
cated hypergeometric functions of the static refractive 
index n± and we have listed them here for completeness 
only. We emphasize that, for the same reasons as be- 
fore, the no-damping limit 711 — > is not permitted in 
the expansion (|69|) . To study the case of no damping we 
must set 7|| = in Eq. (|68p and repeat the asymptotic 
analysis. Then, we get 




a 



Zuj-c 



(70) 



(ju = 1 j 

2n±(n± 



2n 3 ± 



4n 2 ± 



6nj_ 



57ij_ (n± + iy 



Interestingly, the expansions (|66|) and (f69|) derived with 
damping, 7|| 7^ 0, contain fractional powers of Z which 
are not present in the absence of damping, cf. Eqs. (|67p 
and (|70p . This can be traced back to the behaviour of 
a(tu) at lu = 0, which varies depending on the presence 
(or absence) of dissipation in the material. Note that be- 
cause of the presence of the fractional powers of Z in the 
asymptotic expansions ([66)) and ([69]) the significance of 
next-to-leading order terms is much greater there than in 
expansions (|67|) and (f70|) . Because these next-to- leading 
order terms depend on parameters describing losses in 
the material, we conclude that damping plays an impor- 
tant role in the Casimir-Polder interaction between an 
atom and a lossy conductor. This is in contrast to the 
case of an atom interacting with an absorptive insulator 
where even in the presence of damping there are no frac- 
tional powers of Z present in the asymptotic expansion of 
the energy shift in the retarded regime [l[. For lossy di- 
electrics the leading-order behaviour of the energy shift is 
the same as for a non-absorptive dielectric and the next- 
to-leading order term, which depends on the absorption 
parameter, is indeed a small correction. In Fig. [5] we 
illustrate numerically the impact of damping for the case 
of a material that conducts parallel to the surface but 
insulates perpendicular to it. 

Comparing Eqs. ([66]) and ([69]) together with Eqs. ([67]) 
and (|70|) we see that in the case when the response of the 
material in the z direction is taken to mimic a dielec- 
tric rather than a conductor, the quantities pertaining 
to ej_ appear earlier in the asymptotic series. Note that 
if en and e± both describe dielectrics, then already the 
leading-order term of the asymptotic expansion depends 
on both n\\ and nj_, where n_L,|| are the static refractive 
indices i.e. ,1 = e± t \\ (0). Thus, we conclude that in the 
retarded limit the anisotropy of the material is important 
in the case of an atom interacting with a dielectric but 
not so much if it is a conductor. 



3. Case wp|| S> upi 

We now turn to study in greater detail the scenario of 
an atom interacting with a material whose conductivity 
in the direction parallel to the interface far exceeds that 
in the direction perpendicular to the interface. In partic- 
ular, we are interested in the behaviour of the energy shift 
in the non-retarded regime, where the impact of the ma- 
terial's anisotropy is greatest. We assume u)p\\ » cjpj^ 
and thereby effectively set a± =0. We also neglect 
damping, that is, we take 7y = 0. With these simpli- 
fications the shift is still given by the Eq. ([55)) but the 
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FIG. 2. (Color online) The dimensionless F\\ of Eq. (JBDJ that 
enters the shift (|59[) for an electromagnetic response of the 
form (|68p , as a function of the dimensionless parameter Zui mg , 
i.e. the distance from the surface measured in atomic transi- 
tion wavelengths. The solid straight line shows the retarded 
limit for the perfect reflector, which is the leading-order term 
in the asymptotic expansion of F\\, Eqs. ()69[) and (|70|) . In 
the presence of damping (red, dot-dashed) i.e. for non-zero 
711 , even for relatively large values of 2u mg , the corrections 
to the leading-order behaviour are still appreciable, especially 
for poor conductors. Without damping (blue, dashed) they 
are much smaller. 
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(73) 



reflection coefficients are now significantly simpler 



rTE = =5- \ + x 2 - x) - I, 



fTM 



2x 
2xy 



-Pll 
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(71) 
(72) 



where we have abbreviated Wpii = Wpii/w mff . The ?/ inte- 
gral is now elementary and can be calculated analytically, 



where the coefficients F p are understood to enter Eq. 
to give the energy shift. The expressions in Eq. ([73)) 
may look formidable but in fact they are given entirely in 
terms of elementary functions. More importantly, the co- 
efficients F p are given as fast-converging one-dimensional 
integrals [TtJ that can almost effortlessly be computed 
numerically, and the subtleties of the asymptotic anal- 
ysis of two-dimensional integrals need not concern us 

here [lU. Note that terms containing 

in principle become imaginary for wpi 



can 



> 1, although 
in practice this is unlikely to happen. As a typical ex- 
ample let us consider graphite with a plasma frequency 
of the order of uj p w 0.7 x 10 15 Hz [lj| and a rubid- 
ium 87 Rb atom which has its strongest dipole transi- 
tion (5 2 S 1/2 -> 5 2 P 3/2 ) at u mo w 2.4 x 10 15 Hz [H. 



Then, we have wpn = wp/w?, 



0.29. However, even 



if Wpy > 1 the expressions in Eq. (|73[) actually remain 
real, but feature inverse hyperbolic functions through 
arctan(iz) = iarctanh(z). 



Let us now consider the nonretarded limit of Eq. (f73l) . 
In order to obtain the asymptotic expansion we note that 
in the limit Zuj rng — > the exponential in the integrand 
approaches unity and the integral diverges. Thus we re- 
place the integrand by its large- x behaviour and retain 
only the positive powers of x. Then, carrying out the 
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integrations we obtain 



Zuj. 



1 8 2£5p|| 4£iL 



2w P || 



Zuj. 



-arccot 



1 



(74) 



1-2 P|| 



aJpii 2w 



-arccot 



l - uL 



(75) 



Inserted into Eq. ([59]) these results show that to lead- 
ing order the shifts depend on distance as Z~ 3 , as one 
would expect for electrostatic interactions. We note 
that the leading-order contributions each contain a tDpy- 
independent term, which is the only one surviving the 
limit (Dp || — > oo where it yields the perfect reflector re- 
sult, as naively expected. Thus, one could see the results 
([TH [75]) as those for a perfect reflector but amended by 
additional material-dependent corrections that are of the 
same order in the parameter Zuj mg . Note also that both 
and F 1 - vanish in the vacuum limit wpn — > 0, as they 
should. In Fig. Owe compare the asymptotic expansions 
in Eqs. (|74l [75]) with the results of the exact numerical 
integration of Eq. (|73l) . 
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ing part of the integrand with its small- a; behaviour. In- 
tegrating the resulting expression term by term we arrive 
at 
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(76) 



(77) 



This expansion of course coincides with the one in 
Eq. (|70p in the limit n± — >• 1. A comparison of asymp- 
totic and numerical results in the retarded limit is shown 
in Fig. |H 




FIG. 4. (Color online) For large Zuj mg , the coefficients F p 
according to Eq. (|75[) approach constants, which implies that 
the energy shift in Eq. (|59[) behaves as Z~^ in the retarded 
regime. The solid lines represent results of the exact numeri- 
cal integration of Eq. (|73[) whereas the dashed lines have been 
calculated from the corresponding approximate formulae in 
Eqs. (|76l I77[) . Note that the corrections to the leading-order 
terms may be significant for small values of Up/oj mg . 



FIG. 3. (Color online) For small Zui mg we choose to plot 
F p /(Zu) m g) which according to Eq. (|75[) approaches a con- 
stant as Zuj mg 0, so that the energy shift in Eq. ()59[) be- 
haves as Z~ 3 in the nonretarded limit. The solid lines repre- 
sent the results of the exact numerical integration of Eq. (J73J) 
whereas the dashed lines come from the corresponding ap- 
proximate formulae in Eqs. (|74l I75[) . 

In the retarded limit the asymptotic expansion of 
Eq. (|T3"|) is obtained by observing that the integrand is 
again strongly damped in the limit Zuj mg — > oo. Thus, 
we separate out the exponential and replace the remain- 



V. CASIMIR-POLDER SHIFTS NEAR AN 
ANISOTROPIC SLAB 

In this section we are going to derive the energy 
shift for an atom close to a slab of material that is an 
anisotropic conductor or dielectric. We will tailor our re- 
sults to the scenario of an atom interacting with a multi- 
layer of graphene. To do so we take the electromagnetic 
response of the slab to be described by 

e|| ( w ) = l +i 3lM i e ±(u) = n 2 ± , (78) 
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that is, we assume that the slab behaves as a lossy con- 
ductor along its surfaces and as a non-dispersive and 
non-absorptive dielectric in the direction normal to its 
surfaces. Such a system loosely mimics a finitely thick 
wall of stacked graphene sheets where electrons are free 
to move along the sheets but their motion across the 
sheets is severely restricted. 



Here the reflection coefficients are now given by 



R 



1 - ?2 e -2fc^L 



(84) 



with r\ being the single-interface coefficients listed in 
Eq. (|62|) . The wave vectors in the dielectric medium k^ d 
are now 



A. Ground state shifts 

The shift of the atomic ground state due to the pres- 
ence of an anisotropic dielectric/conducting slab may be 
written down instantly by replacing the reflection coeffi- 
cients in Eq. (|56|) with those appropriate for a slab, given 
in Eq. fl47j. We get 



AE„ 



87r 2 e 



E 



oo 



dkk I duj 



UJ 2 +0J 2 ng yjk 2 + U} 2 
|2 , ni,2nTM| 



x{[(fc 2 + w 2 )i?™- W 2 i? TE ] U ig \ 2 + 2k 2 R™U ig \ 2 }, 



(79) 



where R, apart from a phase factor, is the reflection co- 
efficient for a slab; expressed in terms of the appropriate 
variables it reads 



R x =f x 



1 - e 



1 — r\e~ 2k =d L 



(80) 



Here r\ are the same as in Eqs. ([571 158]) and the perpen- 
dicular components of the wave vectors in the dielectric 



medium k^ d are 



fcjj 3 = Je^)u: 2 + k 2 , 



e±(i(j) 



y/ e±(iuj)uj 2 + k 2 . 



K zd ~ ■ LUJ, n 



£|i (ixyui 



my 



-1 



1. 



-TM 



J rng A 



e±(ixyuj mg ) 



[ej_(ixyu) mg ) - 1] y 2 + 1. 



Within the system of units we use in this paper (h = 1 = 
c) the combinations Zuj mg and Loj mg are dimensionless 
quantities, so that the functions F\i± are also dimension- 
less. 



1. Nonretarded limit 

Restoring the missing factors of the speed of light c 
in Eq. (fT9")) . just as we did in Sec. IIV A 11 and taking 
the limit c — > oo, we find that when the atom is close to 
the interface the energy-shift of the ground state may be 
approximated by the formula 



1 



iTT 2 e 



E 



dkk 2 e 



-2kZ 



dw TZ(k, lo) 



(85) 



where we have defined 



Tl{k,uj) 



1 — e ~ 2 ^/ e l\/ e -L kL 



i 



i 



Later we are going to evaluate the energy shift numeri- 
cally, whence we wish to rewrite Eq. (j79p in a form which 
is more suitable for numerical analysis. To this end we 
go to polar coordinates, ui = xoj mg cos <j>, k = xui mg sin <j>, 
and then set co$>(j) = y. After a short calculation we 
obtain 



1 _ y± _ 

4 ^ CO 



| 2 + ^I^J 2 



(81) 



with \iAig\ 2 



l mg\ 



\n v mg \ 2 and 



(Zuj mg ) 4 



F" = / dxx 3 I dy ^ m ?> (R™-y 2 R TE 
1 + arjr V 



-2Zu 



(82) 



tdy ( ZUJm ^ fl o ZTM„-2Z UmaX 



v l-0 2R^e 

1 + x^y z 



(83) 



Here eiij_ = e»±(iu>) are evaluated at imaginary fre- 
quency and are completely general. The non-retarded 
limit Zuj mg — > commutes with the mathematical limit 
pertaining to dielectric-conductor transition, u>r||.x 0, 
as well as with the limit of no damping, 711 x 0- I n 
particular, as long as we remain in the non-retarded or 
van der Waals regime, we may apply formula (j85jl just 
as well with a dielectric response as in Eq. (j65[) as with 
that of Eq. (fTSj). 



2. Retarded limit 

In order to find an asymptotic expansion of the energy 
shift in the retarded regime we take the limit Zoj mg — > 
00 in Eq. (|8ip . However, this limit does not commute 
with the dielectric-to-conductor limit cjt|| ,x ~^ nor with 
the no-damping limit 71 m ~^ 0, so that we must specify 
the electromagnetic response of the slab from the outset 
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of the calculation. Let us restrict ourselves to the case 
when the slab has the optical properties given by Eq. 
(|75| . For large Zuj mg the exponential factor in Eqs. (|82"j) 
and (|83[) strongly damps the integrand so that the major 
contribution to the x integral comes from the vicinity 
of the point x = 0. Taylor-expanding the rest of the 
integrand around this point and then carrying out the 
x and y integrations, wc find that, to leading order, the 
functions FH ,J - entering Eq. ([8Tj) behave as 



1 



1 



i + 

8 l.2L 2 <Tf?(0) 8 WO) 



La\\(0) 
16 
1 



FVjj(O) 
3 



In 
In 



1 



La\\(0) 
La JO) 



8 



Lctm(O) 



Lct||(0) - 1 



l-LVjf(0) 



In 



Lo-m(O) 



(86) 



where <t\\(0) 



1/271 1 is the conductivity of the slab 



along its surfaces evaluated at zero frequency. We note 
that, unlike in the case of an atom interacting with an 
anisotropic half-space, cf. Eq. ([69| . the leading term of 
the asymptotic expansion in the retarded regime does 
depend on the optical properties of the slab (although 
it would independent of them if we had considered a 
lossless conductor). This interesting comparison demon- 
strates that cavity QED calculations based on artificial 
constructs such as a half-space, i.e. an infinitely deep con- 
ductor, can be potentially misleading. 

If, instead of a conductor, we consider a slab made 
of a dielectric material (i.e. a material with a dielectric 
constant of the form Eq. (|17p with non-zero wt) then, 
using the same procedure as described at the beginning 
of this section, we find that the leading-order terms of 
the asymptotic expansions of F"'- 1 in the retarded regime 
read 



F" 



F J 



L 

Z 
L 

Z 



e,,(0) 



2^-5^(0) 



1 1' 

4e±(0) ~ 5 
2 1 

To 



(87) 



where en j_(0) = 1 + 



,±\uj — i t u/pii _i/<^x\\ ± * s ^ ne dielectric con- 
stant evaluated at zero frequency. This is in agreement 
with the common knowledge dating back to the work of 
McLachlan [T(| that it is the static response of the mate- 
rial that matters in the retarded regime of the Casimir- 
Polder interaction. By contrast, the dielectric response 
of an anisotropic conductor in Eq. (|78|) . which was used 
to derive the asymptotic expansion in Eq. (|86p . behaves 
in the static limit as 



On the basis of that one may be tempted to expect that 
in the retarded regime the leading-order behaviour of the 
energy shift near a conducting slab with an e(u>) as in 
Eq. ([75)1 should be that of a perfect reflector, but Eq. (j8"B")) 
has shown otherwise. The static dielectric response in 
Eq. (|88p facilitates perfect reflection for a half-space but 
not for a slab, as one can easily check by taking the zero- 
frequency limits of the appropriate reflection coefficients. 
For atoms interacting with conductors it is the static 
conductivity erii(0) that matters, and not en(O). Since 
<7||(0) depends on the damping constant 711 we conclude 
that damping plays an important role in the retarded 
Casimir-Polder interaction between an atom and a real- 
istic conductor. 

As a final remark we note that the leading-order term 
of the asymptotic series in the retarded regime is indepen- 
dent of n± and therefore our results and conclusions triv- 
ially extend to an isotropic conductor. Just as in the case 
of an atom interacting with a half-space, the response of 
the material in the direction normal to the slab's inter- 
faces is to leading order irrelevant in the retarded regime, 
but it does very much matter in the non-retarded regime, 
as shown by Eq. 



e||(w -> 0) 



00, 



ej.(w ->■ 0) 



B. Numerical examples 

We now proceed to give a few numerical examples 
which illustrate the impact of the anisotropy of the mate- 
rial on the Casimir-Polder interaction between a neutral 
atom in its ground state and an anisotropic dielectric 
slab. For numerical purposes Eq. (f8"Tj) is the most suit- 
able. The numerical integration is straightforward as the 
integrals converge quickly because of the damping pro- 
vided by the exponentials, and it can be carried out with 
standard computer algebra packages. This ceases to be 
the case only when the atom is very close to the sur- 
face but under such circumstances our model would not 
be valid anyway because the main effect of interaction 
would then be due to the overlap of the atomic wave 
function with that of the solid and thus of an entirely 
different nature to what has been investigated here. 

To verify the asymptotic analysis of the energy shift of 
Eq. (fSTj) in the retarded regime, Sec. IV A 21 we plot in 
Fig. [5] the exact value of with and without losses in 
the material. In the case of no damping (blue, dashed) 
we observe that as Zu} mg increases F" approaches a con- 
stant value of 1/2, which indicates that, to leading order, 
the energy shift decays as l/Z 4 with the atom-surface 
separation, with a factor of proportionality that is the 
same as in the case of a perfectly reflecting surface. This 
confirms our remarks made in Sec. IV A 21 regarding the 
lossless slab. When the material is absorbing (red, dot- 
dashed) the function F" still approaches a constant value 
i.e. the energy shift still vanishes as l/Z 4 , but the fac- 
tor of proportionality (red, dotted) depends on the static 
conductivity of the slab along its surfaces and its thick- 
ness through the dimcnsionlcss combination Feri|(0). 
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FIG. 5. (Color online) Exact values of f" of Eq. |g2| with 
e(w) as in Eq. (|78[) . with and without the losses in the mate- 
rial. In the case of no damping (blue, dashed) f 1 " approaches a 
constant value of 1/2 for increasing Zui mg , which is how a per- 
fect reflector would behave in the retarded regime. When the 
material is absorbing (red, dot-dashed) F" still approaches a 
constant value (i.e. the energy shift still vanishes as 1/-Z 4 ) 
but the proportionality factor depends on the dimensionless 
combination La^(0), cf. Eq. (1861) . 
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FIG. 6. (Color online) The function F x of Eq. (83]) 
an isotropic conductor (blue, dashed) and for an e(ui) as in 
Eq. (|78|) . where the conductivity of the slab in the direction 
normal to the surface is suppressed and the slab's optical re- 
sponse in this direction is that of a non-dispersive dielectric 
(red, dot-dashed). As Zui mg increases the difference between 
the energy shifts in the two cases vanishes, which confirms 
our conclusion that the anisotropy of the conductor does not 
matter in the retarded regime. 



Next we would like to illustrate the impact of the ma- 
terial's anisotropy on the Casimir-Polder interaction. For 
this purpose we plot in Fig. [S]the function F ± of Eq. 
for two different cases. In the first case we take the ma- 
terial to be an isotropic conductor (blue, dashed), and 
in the second case we assume that the conductivity of 
the slab in the direction normal to the surface is sup- 



pressed and that the slab's optical response in this di- 
rection is that of a non-dispersive dielectric (red, dot- 
dashed). Fig. [5] illustrates that in the retarded regime 
the anisotropy of the material becomes irrelevant. How- 
ever, Fig. [6] is not convenient for looking at the energy 
shift in the non-retarded regime, where it is much more 
convenient to plot F 1 - j 'Zu! mg because this combination 
approaches a constant value for Zuj rng — > 0. We do just 
that in Fig. which shows that the anisotropy of the con- 
ductor strongly affects the Casimir-Polder interaction in 
the non-retarded regime. 
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FIG. 7. (Color online) The function F ± /Zcu mg of Eq. (JHSJ) for 
an isotropic conductor (blue, dashed) and for an anisotropic 
conductor with an e(oj) as in Eq. (|78[) (red, dot-dashed). For 



small distances Zui mg the anisotropy of the material has a 
considerable impact on the Casimir-Polder interaction. In this 
example, the Casimir-Polder force at Zuj mg ~ 1 is reduced by 
about 25% due to the suppression of the conductivity in the 
direction normal to the slab's surface. 



VI. SUMMARY AND CONCLUSIONS 

The quantum propagator of the displacement field 
describes the propagation of photons in the multipo- 
lar formulation of nonrclativistic QED. Modelling an 
anisotropic dielectric as a set of quantized harmonic oscil- 
lators coupled to a heat bath we have successfully derived 
the quantum photon propagator near an anisotropic slab, 
which may be taken to be either a conductor or a dielec- 
tric, or even one type in one direction and another type 
perpendicular to it. We have used the photon propagator 
to work out the energy-level shifts for a variety of scenar- 
ios, but perhaps most importantly for an atom close to 
a slab that in the direction parallel to its surface be- 
haves as lossy conductor and in the direction normal to 
it as a nondispersive and non- absorptive dielectric. We 
have found that the anisotropy of the material matters 
most in the non-retarded regime and that it affects the 
Casimir-Polder interaction considerably when the atom 
is close to the surface. In the opposite case, when the 
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atom is far away from the surface, the significance of the 
material's anisotropy depends on the optical response of 
the material one deals with. For atoms interacting with 
anisotropic conductors the anisotropy does not matter 
to leading order, but it does matter for anisotropic di- 
electrics. This is in contrast to the role of damping in 
the retarded regime: it does not matter for dielectrics 
but it is important for conductors. Our results for the 
energy shift may be used to estimate the Casimir-Poldcr 
force acting on quantum objects trapped close to mul- 
tilayers of graphene or graphite. Our results are par- 
ticularly important for the case of cold molecules whose 
dispersive interactions with surfaces often fall within the 
non-retarded regime where the anisotropy of the material 
strongly affects the Casimir-Polder force. 



a. Nonretarded limit: In the limit of instantaneous 
interaction, 2|w m i|Z — > 0, Eq. (|Aip becomes 



as: 



i 



E 



x /e\\(\uj m i\)e±(\uj m i\) - 1 



327re Z 3 ^ (|w mi |)ex(|u; mi |) + 1 
x (|mL| 2 + 2|^| 2 ). (A3) 

whose real part gives the residue contributions to the 
energy shift of the excited state \i) to leading order as 



AE* 



1 ^ klld^miDej-d^miDI - 1 

~ 32neoZ3 h l^u(l^™l)^(K^T+ 1 l 2 

x (|mL| 2 + 2|^| 2 ). (A4) 
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Appendix A: Excited-state energy shifts near an 
anisotropic half-space 

The additional energy shift of an excited state \i) (and 
the spontaneous decay rate as its imaginary part) can be 
obtained by using formula (|54|) with the propagator of 
Eq. pTj) . The general formula for the additional energy- 
level shift is given by 



A£ t = E 



'ZOO 

3 ' dn e 2iluJ "" lZK 



Thus, in the nonretarded regime the residue contribu- 
tions depend on distance as Z~ 3 and therefore are of the 
same order as the ground-state shifts. 

b. Retarded limit: In the retarded limit Eq. (|A1[) be- 
comes 



~ ImLI 2 + 2 (1 - « 2 ) r™|^| 2 } 



where \^ nl \ 2 = \[x x mi \ 2 + |/i;y 2 , and 
is the modulus of the transition frequency between the 
states \i) and |m). The dimcnsionless variable k depends 

on the original q» through k = 



u Li - <l\\/\umi\ and its 



contour of integration runs from k = 1 along the real axis 
to k = and then up along the imaginary axis to K = ioo. 
The reflection coefficients expressed as functions of k are 



1 



+ y/ejj -1+K 



2' 



-TM 



(«) 



y?e\\e±n - \/e± - 1 + 



^/^elK + Vel - 1 + ■ 



(A2) 



Here the dielectric permittivities are evaluated at the 
real frequency |w ml |, i.e. e a = e CT (|w rni |). The asymp- 
totic analysis of Eq. (|A1[) can be carried out by exactly 
the same methods as those described in Sec. VI B of 
Ref. [l|. Therefore we merely quote the results. 
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(A5) 



with the parallel refractive index nii(|u;. 



(Al) 



\fe\\ (\uj m i\). Thus, the energy-level shifts, given by the 
real part of the above expression, read 
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m<z 1 II v ' 1 
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(A6) 



Appendix B: Spontaneous decay rates near an 
anisotropic half-space 

In the non-retarded limit the spontaneous decay rates 
are given by the imaginary part of Eq. (|A3[) 
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1 y-v Im[^/ey(|a; m i|)ei.(|a; T ; 



x (|mLI 2 + 2|m^I 2 ), (bi) 
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whereas in the retarded limit by the imaginary part of 
Eq. |(33l: 



Ar rot = _ 



27re \n\\(uj mi ) + 1| 2 



+2 



(|n||(|w ml |)| 2 -l)sin(2|w ml |Z) 
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(2k M |Z) 2 

with the parallel refractive index n||(|w m i|) = 
*/eji (\uj m i\). The behaviour of the excited-state shifts 



and spontaneous decay rates in the retarded regime is 
very much different from the energy shift of the ground 
state in the far zone because they oscillate. This differ- 
ence can be explained by the fact that the ground-state 
shifts are caused by virtual photons whereas an excited 
atom may emit real photons. Therefore, the shift of an 
excited level bears a close analogy to a classical dipole os- 
cillating with the frequency \ui m i\ near an interface. The 
analogy persists for an anisotropic half-space. Note that 
the results (|A6|) and (|B2[) do not, to leading order, de- 
pend on e±(u). This is plausible if one notes that the 
electric field of a radiating dipole that is sufficiently far 
away from a surface is to a good approximation parallel 
to the surface. Thus it is insensitive to the response of 
the material in the direction perpendicular to the surface. 
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